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Abstract 

We develop theoretical formulas for the prediction of the rupture of systems which are known to exhibit a critical 
behavior, based solely on the knowledge of the early time evolution of an observable, such as the acoustic emission rate as 
a function of time or of stress. From the parameterization of such early time evolution in terms of a low-order polynomial, 
we use the functional renormalization approach introduced by Yukalov and Gluzman to transform this polynomial into 
a function which is asymptotically a power law. The value of the critical time tc, conditioned on the assumption that tc 
exists, is thus determined from the knowledge of the coefficients of the polynomials. We test with success this prediction 
scheme with respect to the order of the polynomials and as a function of noise. 



Q 1 Introduction 

O . 

To what extent can the material failure of a mechanical system under stress be forecasted? This question has enormous 
3 technological interest for its economic and human consequences, especially in the automobile, naval, aeronautics and space 
^ industries as well as in the sensitive chemical and nuclear industries to cite a few among many examples. 
1^ If rupture occurs brutally without precursors, prediction is impossible. In contrast, Mogi noticed that, for experiments on 
^ -a variety of materials, the larger the heterogeneity of the material, the stronger and more useful are the precursors to rupture 
Q .[2]. For a long time, the lapanese research effort for earthquake prediction and risk assessment was based on this very idea 

>::i- 

V~j ■ These empirical facts have been put on a firm theoretical basis by using various simplified mechanical models of failure 
of heterogeneous materials which showed that, increasing the disorder changes rupture from first-order (abrupt) to critical 

5^ (continuous with power-law properties) [Q]. By the term "disorder', we refer to heterogeneity in material properties (elastic 
coefficients and rupture thresholds) as well as inhomogeneous pre-stresses. In the presence of long-range elasticity, disorder 
is found to be always relevant leading to a critical rupture. However, the disorder controls the width of the critical region 
[||]. The smaller it is, the smaller will be the critical region, which may become too small to play any role in practice. 

The potential for predicting rupture thus seems associated with its critical nature. Let us first review past works which 
support this concept of critical rupture. Theoretically, the concept that rupture is critical was first formulated by Vere-Iones 
using critical branching models and Allegre et al. [Q] using the percolation model and real-space renormalization group 
(see ref. for a general presentation). The Russian school has also extensively developed this concept |Q p^ ll]. One 



of us and co-workers have introduced a statistical two-dimensional model of dynamically evolving damage [|12|], which 
exhibits the critical time-to-failure dependence of the energy released up to the rupture. Based on a precise numerical 
description of the many growing interacting micro-cracks with a spatio-temporal organization which is a function of the 
stress-dependent damage law, one finds that, under a step-function stress loading, the total rate of damage increases on 



average as a power-law of the time-to-failure. In this model, rupture is a "critical point" in the statistical physics sense |j8|] 
and occurs as the culmination of the progressive nucleation, growth and fusion between microcracks, leading to a fractal 
network of cracks. This simple model has since then been found to describe quantitatively the experiments on the electric 
breakdown of insulator-conducting composites [|l^] and the damage by electromigration of poly crystalline metal films [14]. 
This led to the proposal and the empirical test on real engineering composite structures that failure in fiber composites is 
a genuine "critical" point [15]. This critical behavior may correspond to an acceleration of the rate of energy release or 
to a deceleration, depending on the nature and range of the stress transfer mechanism and on the loading procedure. A 
generalization is to extend the power-law behavior of the time-to-failure analysis by including corrections in the form of 
log -periodic modulations [|l^]. Log -periodicity is the hallmark of a hierarchy of discrete characteristic scales in the rupture 
process. Mathematically, it corresponds to adding an imaginary part to the exponent z (defined below). Intuitively, the 
log-periodic oscillations are oscillations that are periodic in the logarithm of the time-to-failure and thus corresponds to 
an accelerating frequency modulation as the critical time is approached. This acceleration of alternating ups and downs 
accounts for the succession of damage and quiescent phases self-organizing and culminating in the rupture. 

Following this work [15], this method has been used extensively by the French Aerospace company Aerospatiale on 
pressure tanks made of kevlar-matrix and carbon-matrix composites used on the European Ariane 4 and 5 rockets. In this 
application, the method consists in recording acoustic emissions under constant stress rate. The acoustic emission energy as 
a function of stress is fitted by the critical theory mentioned above. One of the parameters is the time of failure and the fit 
thus provides a "prediction" when the sample is not brought to failure in the first test. The results indicate that a precision of 
a few percent in the determination of the stress at rupture is typically obtained using acoustic emission recorded about 20% 
below the stress for rupture. We now have a better understanding of the conditions, the mathematical properties and physical 
mechanisms at the basis of log-periodic structures [[T^. 

The numerical simulations of Sahimi and Arbabi [ p^ ] have confirmed that, near the global failure point, the cumulative 
elastic energy released during fracturing of heterogeneous solids with long-range elastic interactions exhibit a critical be- 
havior with observable log-periodic corrections. Molecular dynamics simulations of the geometry of fracture patterns in a 
dilute elastic network give similar results [jl^] : under a uniform strain which drives the fracture to develop by the growth and 
coalescence of the vacancy clusters in the network, there exists a characteristic time at which a dynamical transition occurs 
with a power law divergence of the average cluster size. The cluster growth near the critical time also exhibits spatial scaling 
in addition to the temporal scaling, namely as fracture develops with time, the connectivity length of the clusters increases 
and diverges at tc- Recent experiments on the rupture of fiber-glass composites have also confirmed the critical scenario [[l9|]. 
Johansen and Sornette [ 20 j have recently re-analyzed the acoustic emissions recorded during the pressurization of spherical 
tanks of kevlar or carbon fibers impregnated in a resin matrix wrapped up around a thin metallic liner (steel or titanium) 
fabricated and instrumented by Aerospatiale-Matra Inc. These experiments were performed as part of a routine industrial 
procedure, which tests the quality of the tanks prior to shipment. It was found that the seven acoustic emission recordings 
of seven pressure tanks which was brought to rupture exhibit clear acceleration in agreement with a power-law "divergence" 
expected from the critical point theory. 

At the same time, it became tempting [21 ] to apply similar considerations to earthquakes. Indeed, over the years there has 
been a growing evidence that a significant proportion of large and great earthquakes are preceded by a period of accelerating 
seismic activity of moderate-sized earthquakes. These moderate earthquakes occur during the years to decades prior to the 
occurrence of the large or great events and over a region much larger than its rupture zone. Sornette and Sammis 

[gl]l iden- 
tified a specific measurable signature of this criticality in terms of a power-law acceleration of the Benioff strain previously 



interpreted as an exponential mechanical-damage rate [E2, 23]. The combined observational and simulation evidence now 



seems to confirm that the period of increased moment release in moderate earthquakes signals the establishment of long- 
wavelength correlations in the regional stress field. Large or great earthquakes appear to dissipate a sufficient proportion 
of the accumulated regional strain to destroy these long wavelength stress correlations [24]. They can thus be considered 
as different from smaller earthquakes. According to this model, large earthquakes are not just scaled-up version of small 
earthquakes but play a special role as "critical points" [25, 26]. Recent extensions to the intermediate scale of rockbursts in 
deep mines confirm the picture [27]. 

To summarize these works, there is a rather strong evidence that rupture in heterogeneous media is critical in the sense 



of statistical physics [||]. To what extent can the critical rupture concept be used to predict rupture? Voight noticed in 
an exciting precursory work that many systems fail by exhibiting a typical relationship relating the second time derivative 



d Q./dt of some observable to a positive power of VL itself [28, M]. He then used this relationship to attempt predictions 



of failures in various materials and of volcanic eruptions. Basically, the relationship he postulates is nothing but a power law 
time-to-failure evolution of the observable Q. 

n{t) = {t, - ty , (1) 

where tc is the critical time of rupture, z is a critical exponent and A a numerical amplitude. By differentiating twice and 
eliminating time, we indeed get d^^/dt^ = B , where B = z{z + 1) 

As we said above, the critical rupture concept establishes theoretically this law (|]) as the result of the cooperative 
organization of precursory damage preparing the global rupture. If an observable such as the rate of acoustic emissions 
radiated during loading exhibit an acceleration close to rupture of the form ([l]) as documented in several experiments [15 
19, 2C], it is clear that one can try to fit the data by (|l]) and get a prediction from the determination of t, 



In practice, the problem is that the fit of a simple power law (|l]) to a noisy data is rather unstable, so much so that often 
no fits can be found [|30t]. The problem comes from the fact that only close to tc (in some relevant time units) can the power 
law be clearly distinguished from other parametric accelerating functional forms, such as an exponential. As an exponential 
has no critical time, this leads to an ill-defined rupture time. The fundamental limitation in using this prediction scheme is 
thus that tc is determined only when data is used up to very close to tc- This is thus far from the prediction goal to infer tc 
from a distance! One possibility to improve the reliability and range of power law fits far from tc has been proposed, based 
on log-periodic corrections to power laws [15|, 1^, 20, 21, |Z/ 



In the present paper, we propose a different approach to the prediction of tc- In the next section, we give the gist of 
the method and formulate the problem in precise terms. We also provide a brief summary of the functional renormalization 
method. In section 3, we test the method in the situation where both the polynomial expansion at early times and the value 
of tc is known, to see how the time-dependence is reconstructed by the functional renormalization. In section 4, we present 
the genuine prediction scheme which determines tc solely from the knowledge of the first few coefficients of the polynomial 
fit at early times and the assumption that the late time dynamics is of the power law form (|l|). In particular, we test how the 
precision of the prediction improves by adding more terms in the polynomial expansion. Section 5 presents numerical tests 
of our method compared to direct fits by a power law for noisy data and section 6 concludes. 



2 Formulation of the problem and outline of the method 

The gist of the method is as follows. We assume that we are able to parameterize the early time evolution of an observable, 
such as the acoustic emission rate as a function of time or of stress by fitting the data to a low-order polynomial Q{t) = 
qq + ait + a2t'^ + .... We then use the sophisticated functional renormalization approach introduced by Yukalov and Yukalov 



and Gluzman [31]-[|42|]. We use the version developed in [38] that transforms a polynomial into an analytic effective sum 



with an asymptotic power law of the form ([I|) close to some tc to be determined. The value of the critical time tc, conditioned 
on the assumption that tc exists, is thus determined from the knowledge of the coefficients oq, ai, ••• 

In order to test our proposed scheme, we compare the results of our method to the exact evolution with time of the 



macroscopic crack length in an exactly solvable model of rupture with damage [|43|, |44|]. In this self-consistent theory, the 
growth of a single macroscopic crack is controlled by cumulative damage dependent on stress history. The damage D 
accumulates according to the equation dD/dt oc cr™, where a is the local stress proportional to the globally applied stress 
(To but which takes into account the distortion due to the crack and m is the damage exponent. The law describing the growth 
of the crack, i.e., the dynamics of its half-length a{t), is obtained from the following self-consistent condition: the time it 
takes from a point at the some distance L from the crack tip at time r for its damage to reach the rupture threshold D* is 
exactly equal to the time taken for the crack to grow from its size at time r by an increment L so that its tip reaches the point 



exactly when it ruptures. For m = I, the full solution is known [43, 44] 



and is indeed of the form (|l|) close to = 3D* /2aQ with z = —1. 

For the convenience of notations, we work in the sequel with dimensionless variables / = o/ao and |^ t ^ t, so that 
the solution becomes 

/exact (i) = 7 • 0) 

cost 

Let us assume that we have access only to the small time dynamics, captured mathematically by the first terms in the 
expansion of the solution @: 

61 

fa{t) = l + ait^ + a2t^ + a3t^ + ... with ai = 1/2, 03 = 5/24 and ag = — . (4) 

Since the cosine function is even in t, only even power of t are present. Knowledge of only the first few terms in (^) is the 
relevant situation for instance in an experiment in which the early acoustic emissions are recorded and one would like to 
infer the subsequent evolution. 



To make this paper self-consistent, we first outline the method we use, which is a direct adaptation of ref. |38]. The 
complete mathematical foundation can be found in earlier publications Assume that we are interested in a function 

(j){x) of a real variable x. Let perturbation theory give for this function the perturbative approximations pk{x) with k = 
0,1,2, enumerating the approximation order. Define the algebraic transform Fk{x, s) = x'^pj^{x). This transform changes 
the powers of the series pk{x), thus changing its convergence properties. As a result, the approximation order effectively 
increases from ktok + s. The inverse transform is Pk{x) = x~''^Fk{x, s). Define the expansion function x = x{f, s) by the 
equation Fq{x, s) = /, where Fq is the first available approximation and / is a new variable. Substituting s) back to 
Fk, we get Ukif, s) = Fk{x{f, s), s). The transformation inverse to the latter reads Fk{x, s) = yj^{FQ{x, s), s). 

Consider the family {y^} as a dynamical system in discrete time k, the order of the approximations. The trajectory 
{Vkifi s)} of this dynamical system is, by construction, bijective to the approximation sequence s)}. This system 

can thus be called the approximation cascade. The next step is to embed the discrete sequence {yk{f, s)} into a continuous 
sequence {y(r, /, s)} with r G [0,+oo]. Thus, the family {y{T,f,s) : r G [0,+oo]} composes a dynamical system 
with continuous time, whose trajectory passes through all points of the approximation cascade trajectory. Such a system 
can thus be called the approximation flow. The evolution equation for a flow can be presented in the functional form 
y{T + T',f,s) = y{T,y{T' , f,s),s). We call this equation the self-similarity relation, which is the central concept of 
our approach. In this framework, the motion occurs in the space of approximations, where self-similarity is a necessary 
condition for convergence as a function of "time" defined as the approximation number. The evolution equation for the 
approximation flow can be rewritten in the differential form and then integrated over time between k and some k* . The 
point k* is to be chosen to provide the best approximation F^j^^{x, s) = y{k* , Fq{x, s), s) for the minimal time k* — k. 
The cascade velocity Vk{y, s) in the vicinity of the time k may be presented by the Euler discretization of the flow velocity 
giving Vkif, s) = Vfc(x(/, s), s), with Vk{x{f, s), s) = Fk^i{x, s) — Fk{x, s). The integral form of the evolution equation 

"'■^^ ^ = k*-k, (5) 

where = F^if, s) and F^_^_^ = F^^^{x, s). The approximation F^j^^ must be reached during the minimal time. When no 
additional constraints are imposed, then the minimal time corresponds, evidently, to one step: k* = k + 1. Finding F^{x, s) 
from (||) and using the inverse transform leads to the self-similar approximation p*^{x, s) = x~'^F^{x, s). 

Now, by means of the substitution s — > s^, we have to introduce control functions s^ which can govern the convergence 
of the sequence Sfc)}- Following the standard procedure [^]-[42] and similar to the steps described above, we may 



construct an approximation cascade with its corresponding cascade velocity. Convergence of an approximation sequence is, 
in the language of dynamical theory, the same as the existence of an attracting fixed point for the corresponding approxima- 
tion cascade. If the cascade trajectory tends to a fixed point, this means that the flow velocity goes to zero as "time" k goes 
to infinity. In practice, we cannot, of course, reach the limit k ^ 00, and have to stop at a finite (approximation order) k. 
Then, the condition to be as close to a fixed point as possible is the minimum of the velocity. 



After the control functions are found from the minimal-velocity condition, we substitute them into and obtain the 
final expression = p^x, Sk{x)) for the self-similar approximation of the sought function. The practical way of using 
the minimal-velocity condition is to express it as a minimal difference condition. To check whether the obtained sequence 
{/^(x)} converges, we have to analyze whether the corresponding mapping is contracting. The mapping related to the 



sequence {/^} is constructed in the standard way [31]-[p2[] and the contraction, or stability, is analyzed by calculating the 
mapping multipliers. Ref. [38] has shown how this method of algebraic self-similar renormalization works to obtain accurate 
estimations of the critical behavior of a large variety of physical systems, starting from virial-type or perturbation expansions 
containing only second-order terms and derived for a region far from the critical point. 

In the sequel, we translate this formalism to the case where x is now a real time, k is the order of the polynomial fit to 
the early time of the signal which plays also the role of "time" for the dynamical flow in the functional space. The goal is 
to describe as accurately as possible the finite-time singularity, which is equivalent to a critical point in the time domain. 
Before investigating the predictive power of our approach, we first investigate the possibility of reconstructing as faithfully 
as possible the time evolution based on the knowledge of the singularity. 



3 Reconstruction of the crack dynamics from the knowledge of the small time dynamics 
and the position of a finite-time singularity 



Here, we assume in addition that we know the position of the singularity att = ^ and its exponent, i.e., that 

1 



m 



tt/2 - t 



for t^TT/2. 



(6) 



Note that the amplitude A of the pole is exactly equal to 1. We apply the Yukalov-Gluzman technique in its version related 
to crossover phenomena [|39|]-[p2|] , to obtain the best approximation of f{t) for arbitrary times based only on the knowledge 
at small times and @ at times close to rupture. 

3.1 Using only fait) = 1 + ait^ and (@) 

From the expansion fa{t) = 1 + ait'^, the Yukalov-Gluzman method allows us to build the approximant 

-0 



flit) 



(exp ( 



air 



(7) 



with two unknown parameters b and (3, to be determined solely by demanding the existence of a crossover. The exponent 
/3 is determined from the condition that fait) = 1 + ait^ must cross-over to (Q). This gives /? = 1 for the simple pole (^. 
The coefficient b is obtained from the condition of a pole at the known critical point tc = vr/2, which reads l//i(ic) = 0. 
Solving for b, we finally get the first-order approximant 



flit) 



exp [-ait^^ - ^ exp {-aitfj 



-1 



(8) 



This approximant gives an amplitude Ai = 0.834 for the simple pole, only 17% off the exact value A = \. Figure [I| 
plots in logarithmic scale the relative errors between the approximant /i and the exact expression (|3|) as a function of time t. 

3.2 Using fait) = 1 + ait^ + aat^ and ® 

A similar procedure as in the previous case gives the second-order approximant 



fm 



exp [ —ait exp 



02^2 

ai 
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exp —ait^ exp 
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n -1 



(9) 



This second-order approximant can be improved by introducing an additional control parameter r (optimal effective 
time k* — k) such that we can enforce the condition that the coefficient 02 is preserved in the renormalization procedure. In 
other words, the expansion of /2(t) is now imposed to have the same coefficient of its power t^. The expression @ is thus 
modified into 



/2W 



exp 



-a\t exp — T t 



+2 / '^2 ,2 
-ait^ exp ( — T 

ai 



(10) 



The amplitude of the simple pole predicted by this approximant is found equal to 0.898, only 10% off from the exact 
value 1. The condition that the coefficient of the power in the expansion of (10) is equal to 02 = 5/24 given by 
to the improved second approximant of the form (|I0|) with 



leads 



1 



2a2 



(11) 



Figure [I] plots in logarithmic scale the relative errors between the approximant /2 and the exact expression (||) as a 
function of time t. 

3.3 Using fait) = 1 + ait^ + 02 + a^t^ and (g) 

An extension one step further in the Yukalov-Gluzman procedure gives 



, .2 I ^2 ,2 ( 03 ,2 

exp ait exp — t ti exp [ — t T2 



ai 



a2 



-13 



(12) 



where the two control parameters ri and T2 are to be determined from the condition that 02 and 03 are conserved by the 
renormaUzation. We find 



1 



2a2 



and 



1 



T2 



"2^1 



+ aia2Ti + 



«3n \ 2ai 

The coefficient ag is found from the condition on the critical point: 
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(13) 
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(15) 



The solution then reads 
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(16) 



This approximant gives an amplitude ^3 = 0.9661 for the simple pole, only 3.4% off the exact value A = 1. The quality of 
the reconstruction of the full function (^ can also be checked by comparing the exact value of the next term (277/8064) = 
0.03435 in the expansion (|) of the function (|). We find 0.034497 t^, corresponding to an error of 0.4%. 

Figure 1 plots in logarithmic scale the relative errors between the approximants /f , and /| and the exact expression 



4 Prediction of the critical time from the knowledge of the small time dynamics 



We now assume the knowledge of the expansion (Q) up to some order, representing for instance the experimental recording 
of an acoustic emission signal up to some stress level. In addition, we assume only the existence of a singularity of the 
form l/{tc — t)^ at some value tc, using the insight from the theory of critical rupture, but do not know a priori neither the 
position of t^,, nor the value of the exponent K. In other words, we assume that we know that the system is bound to break 
but we do not know when and how. Our goal is to attempt to determine the critical time and the functional form of the signal 
on the approach to the critical rupture time from the recording of the early signals. 



4.1 Standard approach from Yukalov and Gluzman [ ]38| ] 

Consider an expansion of some function (p in powers of some variable u given by 



Pk{u) = ^bku'' 

k=0 



with 60 = 1 



(17) 



The method of algebraic self-similar renormalization pq]-[p8[] gives the following general recurrence formula for the ap- 
proximant of order as a function of the expansion Pk-i{u) up to order k — I: 



4>k{u) =Pk-i{u) 



-1 khk k k/s , -s 
1 —u PkLiiu) 



~s/k 



-kjs , X k bk k 
Pk-1 (^) — ^ 



-s/k 



(18) 



where, generally speaking, s = Sk{u), depends on the approximation number and the variable u. 

First, let us estimate the position of the critical point. Using only fi{t) = 1 + ait^, we can write it as the inverse of a 
function that is requested to vanish in order to obtain the singularity. Expanding in powers of t up to first order in t, we get 
the estimation tci = = 1.414, which should be compared with the exact value 7r/2 = 1.571. Including the next order 
from the expansion of cos(t)^^ leads to the improved estimate tc2 = 1-59. 

Let us now obtain the expansion which will be used as a raw material for renormalization. Including the next order in 
f2{t) = flit) + ci2t^, inverting it, and expanding in powers of t up to t^— terms, we obtain an expansion P2{t) for I//2 : 



P2it) 



2 

Y^^kt 

k=0 



2k 



(19) 



1, 61 = -1/2, 



1/24, 



Note that, in the initial series, all coefficients in the expansion of cos~^ are positive, giving the worst possible case for 
resummation. In contrast, the coefficients of the inverted series have alternating signs, which may be better for resummation 
[^]. Hereafter, we apply the resummation procedure to the function F = 1// inverse of /. Correspondingly, /* = (l/F)* . 
In order to determine the exponent K, we follow Yukalov and Gluzman and construct the two approximants avail- 



able from the knowledge of the two coefficients oi and 02. They can be readily obtained from the general formula (18), with 
u = t^. The first order approximant is 

l-^i^ 

si 



(20) 



Representing P2{t) as P2{t) = 1 + hit^{l + ^ t^), and applying the general formula to the expression in brackets, we obtain 
the second order approximant 



F^{t) = l + bit'[l 



blS2 



-S2 



Assume further that 



Sl = S2 



(21) 



(22) 



where s is a single control parameter, i.e. limit of the total control function in the critical point, which will play the role of 



the critical index K. As it was explained in []38|], such an assumption is well justified in the vicinity of a stable fixed point. 
We impose the condition of the existence of a critical point, which delivers two equations for and s: 



F^itc,s) = and F;{tc,s) 







(23) 



The condition of maximum stability of the renormalization amounts to imposing that the difference Fg* — be a minimum 
with respect to the set of parameters. The minimization of the difference is automatically satisfied when (^) holds, since 
the difference reaches it smallest possible value, namely zero. 

The vanishing of F^ given by (20) gives = s/ai. The second condition F2 = with (^TJ) yields the estimation 
s = 1.258 for the critical index, only 26% off the true value equal to 1. The critical time is given numerically by tc = 1.586, 
very close to the exact value 7r/2 = 1.5708. Note that, as is often found in critical phenomena, an error of less than 1% in the 
location of the critical point is associated with a much larger error of about 26% on the exponent. In the scheme presented 
above, it was possible to find s and tc from to separate equations. 

Such convenience does not hold for another approximation scheme, presented below, which has however other advan- 
tages such as simplicity. In order to separate the variables, we need an initial guess either for s or for tc- Such an initial guess 
is provided naturally by the analog of a mean-field approximation. 



4.2 Alternative approach: expansion around a "mean-field" approximation 
4.2.1 Expansion to the same order as above 

An alternative and more transparent approach is first to minimize the distance between approximants and then to verify 



s/bi, i.e., s = bit^. 



that (g3|) holds. This is the reverse order to the previous scheme that solves ( |23| ) which then automatically ensures that the 
distance between the two successive approximant is minimized. 

In practice, this is implemented as follows. From (|^), we see that the critical point is located at 
Using only 

pi{t) = l + ht\ (24) 

we estimate t"^ = —1/bi, which then yields s = —1. Note that this value —1 for the exponent always holds for any value of 
hi. This exponent thus plays a role analog to a mean-field approximation in statistical physics. The fact that the exponent is 
in the present case equal to the exact value is a mere coincidence. 
In the next order, 

P2{t) = I + bit^ + b2t^ . 

The two approximants can be derived directly from the general formulas: 



(25) 
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-S2/2 



We assume in addition that 



S2_ 

2 



(26) 



(27) 



(28) 



thus eliminating a trivial dependence of the control parameter s on the approximation number k. Note that (^ ) is different 
from ( |22| ) because we use a different sequence of approximations /i, /2, .... In particular, the condition ([28[) ensures that 
Fi{t) and F2{t) have the same exponent/control parameter. 

In constrast with the previous method of section 4.1, we first minimize the difference F2{t, s) — Fi{t, s), and then verify 
that (23) holds. The difference calculated at the "mean-field" threshold = —1/bi gives 



l + s 



bis 



(29) 



which is exactly zero at 



hi 



-1.167 . 



(30) 



The position of the critical point tc can be re-calculated from the condition F*{tc, s) = 0, which has a non-trivial solution at 



tc=J- 



S 



(31) 



Expanding in powers of we estimate tc ~ 2^ 



1.532. 



Thus, the renormalization scheme used to calculate tc and K corresponds to an expansion around the mean field value 

tl = l/hi and K = 1 in inverse powers of the dimensionless "Froude" number ^ (see [ p6| ] for a definition and use of the 
Froude number in this context of functional renormalization). In summary, we get the predictions tc = 1.528 closer to the 
exact value 7r/2 = 1.5708, and the critical exponent K = 1.167. 

Higher-order expansion Including the next order in /3(t) = /2(i) + CL^t^, inverting it, and expanding in powers of t up 
to t^— terms, we obtain an expansion pz{t) for I//3: 



k=0 



2k 



with 



61 = -1/2, 62 = 1/24, 63 = -1/720 . 
The two higher-order approximants can be written as follows: 



(1 + ht^) 



-2/s, _ 252^41 
S2 



Assume that 



l + bit^ + h2t^ 
S2 S3 



-3/s3_ 363^e]-«3/3 

S3 



s . 



2 3 

The difference F^{t, s) — F2{t, s) calculated at the "mean-field" critical point tg = — l/^i gives 



D2{s) 



-1/s 



b{J s 



1 



bfs 



(32) 
(33) 

(34) 
(35) 
(36) 

(37) 



which has a zero at s = —1.024. The position of the critical point tc can be re-calculated from the condition l//|(ic, s) = 0, 
which has a non-trivial solution at tc = 1.5722 determined from the equation 



{l + b,tl + b2ttf^'-^t 







(s = -1.024) 



(38) 



which wins over the "trivial solution" of 1 + bit'^ + 62*^ = at 1.592. Note that at order 2, from the condition F2*(^c; s) = 0, 
we could only find the trivial solution = 1/bi. It is only at higher order, starting with the order 3 discussed here, that we 
get corrections to the "mean-field" approximation. 



To test the validity of the expansion and the strength of the corrections to the mean-field approximation, let us represent 
s as s = —l + X, substitute it into equation D2{s) = 0, and expand in powers ofX, thus assuming that X is small compared 
to —1. Keeping only terms linear in X, X can be expressed as follows: 



X 



bib2 1 - In 



(I)' 



(39) 



Such an expansion is justified only when <^ 1 and 1 — In is significantly different from zero. In our case, 

^ = 0.067, 1 - In (||) = 2.792 and thus X = -0.024 is small as expected. 

Similarly, let us represent tf. as tf. = —l/bi + C, substitute it into the equation -F|(tc> s) = 0, and expand in powers of 
C, thus assuming that C is small compared to —1/bi. Keeping only the terms linear in C, we get C = 0.458, which leads 
to the value of the critical point tc = 1.568 (close to 1.5722). The expression for C has a simpler structure at the mean-field 
point s = —1: 

1 / bs ^ 1 

bi \bi 62 



C 



1 



61 b2 



b2 



(40) 



with an estimate for the critical time given by 




bj 

'bi 62 62 



b3 



1.563 . 



(41) 



Thus, the renormalization scheme used to calculate and K corresponds to an expansion around the mean field value 



62 



tc = l/bi and K = 1, in powers of two dimensionless parameters and 

In summary, we get the predictions tc = 1.5722 very close to the exact value tt/2 = 1.5708 and K = 1.024 only 2% 
off the exact value 1. 

We conclude that the applicability of the scheme presented in this section crucially depends on existence of typical small 



parameters, while the original scheme of [ p8[ ] does not have such a dependence and can be applied even if critical indices 
deviate strongly from the mean-field values. 



5 Synthetic tests in the presence of noise 

Figure ^ shows one realization of a synthetic noisy data obtained from (||) with a multiplicative noise of variance 10-2. This 
noisy data simulates an experiment recording a signal as a function of time or of increasing strain or stress. Our goal is to 
use this noisy time series up to a maximum value away from tc to guess using the functional renormahzation method what 
is the critical value tc of divergence (theoretically equal to tt/2 = 1.57...). 

Figure ^ shows the inverse of the function in figure ^ as well as four other realizations, which are used in the functional 
renormalization scheme developed in the previous sections to predict tc- The five different symbols shown in figure ^ 
correspond each to a single time series for a specific noise reahzation. The spread around the theoretical cos t formula gives 
a sense of the amplitude of the multiplicative noise. 

Figure ^ compares the prediction skill of our procedure described in section 4.2 to that from a direct fit with a power 
law. Specifically, we plot the predicted value for the critical time tc as a function of the distance to tc, obtained with the two 
schemes. In this goal, we generated 1 ,000 synthetic data sets using expression (^) and modifying it with a multiplicative 
noise of variance 10^'^ corresponding to a standard deviation of 3.3%, i.e., with the equation fi{t) = ^^[1 + where r] 
is a Gaussian white noise with variance lO^'^. To generate curve a) in figure ^ each data set was fitted by the power law 
equation A{tc — t)"^, with (3 = 1 fixed and A and tc determined from the fit in the time interval [0.5; Tiastpoint]- The thick 
line is the average tc taken over the 1 ,000 realizations and the two thin lines are the it one standard deviations. The curve b) 
in figure ^ is the predicted tc obtained from our resummation technique given in section 4.2 which assumes that the exponent 



P is close to 1, using the coefficients bi and 62 of the fit with the expansion 1 + bit"^ + 62^^ to the inverse of each of the 
synthetic data set fi{t). The thick line is the average predicted tc and the two thin lines are the it one standard deviations. 
The horizontal line at tc = vr/2 is the exact theoretical value for the critical time. Up to very close to the critical time, our 
resummation method is clearly superior to the power law fit, even when knowing a priori the value of the exponent. 

Figure |5| compares the prediction skill of the Yukalov-Gluzman method used in section 4. 1 to that of a direct power law 
fit. Again, 1,000 synthetic time series were generated with multiplicative noise of variance 10^'^. The curve a) of figure^ is 
obtained by using the general resummation method of section 4. 1 which does not assume any specific value of the exponent 
f3. We first fit the inverse of each synthetic data set with a parabolic expression 1 + bit"^ + 62*^ and use these coefficients 61 
and &2 to obtain our prediction The two thin lines are the it one standard deviations. The curve b) in figure § is obtained 
by fitting each data set by the power law equation A{tc — where A, [3 and are all three free parameters determined 
from the fit in the time interval [0.5; Tiastpoint]- The thick line is the average tc taken over the 1,000 realizations and the two 
thin lines are the it one standard deviations. The horizontal line at tc = vr/2 is the exact theoretical value for the critical 
time. Note the striking superiority of our resummation method over a direct power law fit. 

Figure |6| is the same as figure ^ but with a larger multiplicative noise of variance 10~^. Figure ^ is the same as figure ^ 
but with a larger multiplicative noise of variance 10~^. 



6 Concluding remarks 

We have tested two methods for the prediction of the critical time of a singular power law behavior and tested their prediction 
skills against the direct determination by a power law fit. Our analysis and the numerical tests convincingly demonstrate the 
value of our approach which provides significant improved forecasting skills. Our tests have however been restricted to a an 
important from a physical viewpoint but still special case of a function which admits an expansion with only even powers of 
t. The most general situation contains also odd powers of t which complicates the situation. We intend to report progress in 
this general case in a forthcoming publication. 

Finally, we wish to comment upon a conceptual understanding of the role of the exponent s used as control functions in 



the functional renormalization schemes used here. Yukalov and Gluzman [ |38| ] first noticed that their functional renormal- 
ization frameworks allowed them to propose a novel physical understanding of critical exponents as being directly related 
to limits of control functions at the critical point. In other words, they appear as physical analogs of the rather abstract 
mathematical objects given by the control functions. A scale invariant formulation using logarithmic variables allows us to 
understand this rather surprising observation: critical exponents can be seen to be determined by the initial conditions of an 
operator of the group of the symmetry of scale invariance [^]. Indeed, a power law function 0(r) of the distance r = tc — t 
to the critical point, has the property of scale invariance which reads 

0(Ar) = A"(/>(r) , (42) 

where A is an arbitrary magnification factor and a is a critical exponent. This equation means that the field i?i)(r) is invari- 
ant under the homothetical transformation r — > Ar; (p A~"(/>. Expression (^) can be transformed by picking up an 
arbitrary reference field 00 and an arbitrary reference parameter tq and introducing the log- variables U = ln{cl) / (po) , s = 



r/ro, U{0) = In (0(ro)/0o) , U{s) — U{0) = as. In these variables, equation ( |42| ) can also be written: 



Uis) = U{s + fi)+g{f,) , (43) 
9{f^) = U{0) - C/(/u) , for any /i . 

Reciprocally, any regular function obeying (^) is necessarily of the shape U (s) — U{0) = as, where 

a = U{1) - U{0) (44) 



is an arbitrary parameter selected by the initial conditions! 
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Figure 1: Relative error [/cxact(i) - /*(*)] //exact (i), for the three approximations formulas (|) for (t), ([l^) for /^t), 
and @ for /Kt), where /exact (i) is given by (|). 
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Figure 2: Noisy data obtained from @ with a multiplicative noise of variance 10^^. The goal is to use this noisy time 
series up to a maximum value away from tc to guess what is the critical value tc of divergence (theoretically equal to 
7r/2 = 1.57...). 
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Figure 3: 1 over the function shown in figure 2 as well as four other realizations, which are used in the functional renormal- 
ization scheme developed here to predict tc- The five different symbols shown in figure || correspond each to a single time 
series for a specific noise reahzation. The spread around the theoretical cos t formula gives a sense of the amplitude of the 
multiplicative noise. 
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Figure 4: This figure shows the predicted value for the critical time tc as a function of the distance to tc, obtained with 
different schemes, in order to evaluate the value of the resummation procedures proposed here. In this goal, 1 ,000 synthetic 
data sets were generated using expression (js]) and modifying it with a multiplicative noise of variance 10~^ corresponding 
to a standard deviation of 3.3%, i.e., with the equation fi{t) = ^^^[1 + r]], where r/ is a Gaussian white noise with variance 
10~^. To generate curve a), each data set was fitted by the power law equation A{tc — t)~l^ , with /3 = 1 fixed and A and tc 
determined from the fit in the time interval [0.5; Tiastpoint]- The thick line is the average taken over the 1,000 realizations 
and the two thin lines are the it one standard deviations. The curve b) is the predicted obtained from our resummation 
technique given in section 4.2 which assumes that the exponent {3 is close to 1, using the coefficients hi and 62 of the fit with 
the expansion 1 + hit^ + ^2*^ to the inverse of each of the synthetic data set fi{t). The thick line is the average predicted 
tc and the two thin lines are the it one standard deviations. The horizontal line at tc = 7r/2 is the exact theoretical value for 
the critical time. Up to very close to the critical time, our resummation method is clearly superior to the power law fit, even 
when knowing a priori the value of the exponent. 
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Figure 5: Same as figure 4 with 1,000 synthetic time series with multiplicative noise of variance 10~^, and with different 
prediction schemes. The curve a) is obtained by using the general resummation method of section 4.1 which does not 
assume any specific value of the exponent /3. We first fit the inverse of each synthetic data set with a parabohc expression 
1 + + 62*^ and use these coefficients 61 and 62 to obtain our prediction tc- The two thin lines are the dr one standard 
deviations. The curve b) is obtained by fitting each data set by the power law equation A{tc — where A, /3 and tc are all 
three free parameters determined from the fit in the time interval [0.5; Tiastpoiut]- The thick line is the average tc taken over 
the 1,000 realizations and the two thin lines are the ± one standard deviations. The horizontal line at = 7r/2 is the exact 
theoretical value for the critical time. Note the striking superiority of our resummation method over a direct power law fit. 




Figure 6: Same as figure 4 but with a multiplicative noise of variance 10 




Figure 7: Same as figure 5 but with a multiplicative noise of variance 10 ^. 



